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ABSTRACT 



The general 



and the bandwidth response 



of the Kaijo Denki PAT 313-1 sonic anemometer/ thermometer 



were examined. Measurements were compared to those of a 
hot wire anemometer system and a one mil thermocouple sys- 
tem. Turbulence data were obtained by all three systems 
simultaneously in two experiments conducted at a beach 
site. The sonic anemome t er / t he rmome t er system proved to 
be a reliable, rugged, easily installed system with a band 
width response of from 0.0 Hz to approximately 12.0 Hz. 
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INTRODUCTION 



I . 



A. PROBLEM 

Turbulent processes in the atmosphere are of interest 
to meteorologists since it is by these processes that effi 
cient vertical and horizontal transfers of energy, heat 
and momentum occur. These transfer processes have been 
experimentally observed to occur within the frequency band 
from 0.01 Hz to 10.0 Hz. An important method for describ- 
ing the role of turbulence is through use of a system of 
instruments capable of measuring over this band. The ana- 
log data obtained is then read into computer systems to be 
digitized into discrete samples and then to be analyzed by 
spectral or other techniques, 

B . OBJECTIVES 

The primary objective of this thesis was to utilize a 
sonic anemometer/thermometer to measure turbulence and 
evaluate the band over which it is capable of measuring. 

A secondary objective was to utilize a Hybrid computer sys 
tern to digitize the analog data and then to use the IBM- 
360/67 computer system for spectral analysis of the digi- 
tized data. A third objective was to compare a sonic 
anemome t e r / the rmome t e r with a hot wire anemometer system 
and a one mil thermocouple system. 



11 



II . GENERAL DESCRIPTION AND OPERATING PROCEDURES 

OF THE 

SONIC ANEMQMETER/THERMOMETER MODEL PAT-313-1 

A, GENERAL 

The sonic anemome t er / thermome t er^‘ measures the three 
velocity components of the wind by detecting accurately 
the transit time difference (transit time sum for tempera- 
ture) between two acoustic pulses propagated in opposite 
directions across a fixed path. The three velocity com- 
ponents of the wind are the longitudinal, lateral and ver- 
tical components denoted by u, v, and w respectively along 
axes X, Y and Z (Figure 1). The w component is obtained 
directly from the vertical axis, but u and v have to be 
computed from the velocity components and measured 

along the horizontal sonic paths A and B which are 120 
degrees apart. 

B. COMPONENT PARTS OF THE SONIC ANEMOMETER/THERMOMETER 
^ • Power Supply /Indicator Unit (PR-1) 

The power supply /indicator unit and its component 
breakdown is pictured in Figures 2, 3, 4. The instrument 

operates on 100/115 volts with + 10% variation at 50 to 
60 Hz. Total power required is 100 watts. All switches 
and controls are located in this unit. Calibration and 
adjustment procedures (discussed later) are performed at 

^Model PAT-313-1 manufactured by the Kaijo Denki Company 
Limited of Tokyo, Japan. 



12 









4 ' 




M 



this unit. Visua3_ meter output (analog form) is available 
for wind measurement on three range scales of 1) 0 to + 1 

m/sec, 2) 0 to + 3 m/sec and 3) 0 to + 10 m/sec (Fig- 

ure 3). Temperature meter output is + IOC from mean air 
temperature which is dialed in on 0 dial (Figure 4). Two 
BNC output jacks are available for tape recording or sig- 
nal monitoring purposes on both wind and temperature units. 

2 . Probe (TR-31) 

The probe has to be employed within the wind direc- 
tion limits indicated in Figure 5. The wind velocity would 
be less than actual if the probe were positioned so that 
the wind was parallel to the path travelled by the sound 
pulses or if the supporting structure of the probe inter- 
rupted the wind flow. 

The probe heads are ceramic type transducers from 
which sound waves are trasmitted at a frequency of lOOKKs. 
The transducers have high directivity and damping to mini- 
mize interference with other channels and are small physi- 
cally in order to avoid blocking the air flow (Figure 6). 
The transducers screw into special holders designed to 
minimize acoustic coupling through the array. All probe 
heads are interchangeable with only minor adjustments re- 
quired in the wind / t emp era t ure units affected by the 
change • 

The probe heads, holders and entire array structure 
have been designed to reduce interference to the wind flow. 
All wiring and connections have been inserted inside the 
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A,B axes along sonic paths. 
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Figure 1. Wind and Sonic Axes. 
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Figure 2. Power Supply/Indicator Unit (PR-1). 
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Figure 3. Wind Unit (PD-11) . 
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Figure 4. Temperature Unit (PS-11). 
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probe tubing for protection from the elements as well as 
for reducing air flow interference (Figure 7). 

3 . Junction Box (QA-33H) 

The junction box contains three sets of pulse os- 
cillators (which generate the lOOKHz signal sent to the 
transducers) and six sets of head amplifiers (which amplify 
by 1000 the lOOKHz signal sent to the transducers). The 
junction box is a watertight unit with connectors at both 
ends (Figure 8) . One hundred meters of cable link the 
junction box and supply /indicator unit and three meters of 
cable link the junction box and probe (Figure 9). The 
junction box must be connected to the probe and power sup- 
ply prior to turning the sonic anemometer on or the output 
transformer will be shorted. 

C. ADVANTAGES/DISADVANTAGES OF THE SONIC ANEMOMETER/ 

THERMOMETER 

The sonic anemome t e r / t hermome t e r has many advantages 
over previously developed wind and temperature measuring 
systems. Recent results by McBean (1972) revealed the 
advantages the sonic anemometer has over a Gill propeller 
anemometer especially with respect to sensor response and 
angle of ^attack (cosine response) in eddy correlation 
measurements. A report by Kaimal aj^ (1971), in which 
they compared the pulse and continuous-wave type sonic 
anemometers, lists many advantages of the sonic anemo- 
meter system over other velocity measuring systems. 
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(3) Top view 
A- come orient 



fb) Side view 




Figure 5. Wind Direction Limits of Probe Type TR-31 
Kaijo Denki (1969) 




Figure 6. Probe Head, Kaijo Denki (1969) 
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Figure 7. Probe Type TR-31. Kaijo Denki (1969) 
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The PAT-313-1 system has no mechanical parts which 
eliminates lag due to mechanical inertia and thus allows 
the sonic to measure very light and sudden changes of the 
wind. The sonic also has high sensitivity and a broad 
range in all three wind components and in the temperature 
component . 

The resolution of the sonic anemometer is + 2 cm/sec 
in a span of 0-10 m/sec and the resolution of the sonic 
thermometer is + O.IC in a span of + IOC centered around 
a mean temperature range of -IOC to 40C, The response is 
linear and so is its sensitivity, (Therefore, the res- 
ponse is the same at 5 cm/sec as it is at 10 cm/sec.) 

The sonic anemometer has a time response in the milli- 
second range as it detects changes in pulse time between 
two sets of transducers. It takes measurements 440 times 
/sec for the wind components and 220 tiraes/sec for the 
temperature component which allows it to detect changes of 
approximately 100 Hz at an accuracy of three percent 
(Kiajo Denki, 1969), 

The sonic anemometer is not affected by exposure to 
the direct rays of the sun or to the rain as would a thermo- 
meter or a hot wire anemometer. It is, however, affected 
by immersion in water, A high moisture, high heat environ- 
ment affects measurement procedures. These have to be 
corrected in order to get true readings because the measur- 
ing circuits are set for an atmosphere of zero percent 
relative humidity at a temperature of 20C. The probe. 
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junction box and connecting cables are waterproofed allow- 
ing for all weather operation of this system. There are 
no easily broken, sensitive elements (such as a hot wire) 
exposed to the weather. It is a rugged system that can 
accept some rough handling. Cable connection points, es- 
pecially where they enter the probe, are sensitive areas. 
Too much twisting of the cables can result in broken wires 
in the connectors. Care must be exercised so as not to 
bend the probe thus changing the length of transmission of 
the acoustic path and causing misalignment of the probe 

heads. Both problems would create serious errors in velo- 

2 

city and temperature measurement. Each of the probe heads 

has a rubber cap over it to protect the ceramic transducer 

element but abrasive handling will rip it. 

The probe has a small bubble level in the top plate 

which is used to insure that it is properly leveled when 

installed. Errors are introduced if the probe is not hori- 
3 

zontal. The probe should be mounted in such a way that 
it is free of vibrations. It measures relative speed of 
the environment with respect to the probe so any movement 
of the probe will introduce errors into the signals. Any 
object placed within two meters of the probe will reflect 
erroneous signals into the transducers. These are recog- 
nizable as noise spikes. 

2 

Kaimal e t al (1971) has an excellent discussion of these 
errors . 

3 

Kaimal e t al (1971) discusses this error. 
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D. MEASUREMENT PROCESS 
1 . Wind Measurement 

As described previously, the sonic anemometer mea- 
sures wind velocity by detecting the transit time differ- 
ence between two acoustic pulses propagated in opposite 
directions across a fixed path. The nomenclature used to 
describe the measuring process for a typical pair of 
transducers is shown in Figure 10. 

The propagation times t^) for the pulses from 

TRj^ to R^ and TR^ to are 



£ 



( 1 ) 



C cos a + V 



d 



and 



t 



2 C cos a - V 



( 2 ) 



d 



The travel time of the sound with the wind is less 

than the travel time (^ 2 ) with the sound against the wind. 
The difference At = t 2 ~ is given by 



At = ^ (3) 

(C cos a) - V/ 
d 



21V, 

d 



(3) 



Angle a is assumed to be small; therefore, cos a 



2 2 

~ 1. One further assumes that C > > V and therefore, 
2 2 

C > > . These approximations result in 




(A) 



which, when solved for gives 




( 5 ) 



2A 



£ 




V 7 ind Velocity 
and 

Sound Transmitting 




KEY ; £ = distance between transmitter and receiver 

C = sound velocity 
V = wind velocity 

V = component of wind velocity perpendicular 
to sound transmitting path 

= component of wind velocity parallel to 
sound transmitting path 

t- = time for sound to travel from to R. 

(with wind) 

t2 = time for sound to travel from T2 to R2 
(against wind) 

TR^, TR2 = transmitting transducers 

Rf , R2 = receiving transducers 

a = angle between transmitting path and dis- 
tressed sound path 

Figure 10 . Schematic of Typical Measuring Process and 
Nomenclature . 
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The wind velocity parallel to the sound transmitting path 

(V^) is therefore determined and is the quantity indicated 

on the metered output of the sonic anemometer as U or V . 

A B 

The values of u and v (the longitudinal and lateral wind 
components of the X and Y axes of Figure 1) have to be com- 
puted from the velocity components U and V • The w com- 

A B 

ponent is taken directly from the output* 

The equations for u and v are: 



u = (cos 3 ~ ~ sin 3) U. + (cos 3 + ~ sin 3) (6) 

/3 A B 



V =-( — cos 3 + sin 3) U.+ ( — cos 3 - sin 3)V 
/3 ^ /3 ® 



(7) 



where 



3 = tan 



-1 ''b - "a 



/3 (Vj + l.^) 



Overbar denotes time average over a suitably long 
period of time Kaimal e£ aj^ (1971)* 

2 • Temperature Measurement 

Temperature is measured in a manner similar to the 
wind measurement except that the propagation times 
are added instead of subtracted. The w component trans- 
ducers provide the data used for this measurement. The de- 



fining expression for temperature is 

25,(C^ - V 

n 



t = t^+ » 



2 2 



( 8 ) 
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2 2 2 2 
Again assuming C ^ ^ ^ ^ ^ ^n 

equation (8) reduces to: 



t = 



1 % 



( 9 ) 



The speed of sound in air is expressed as: 

C = 20.06 • [T + (1 + 0.32 (10) 

where T is in degrees absolute (K) , e is the vapor pressure 
in atmospheres and p is the atmospheric pressure in atmos- 
pheres. From (9) and (10) 



X = 1 (11) 

402 .68(1 + 0.32 -) 

P 

Furthermore, set T„„ = T(1 + 0.32 — ) where T„., is 

b V P b V 

sound virtual temperature (taken as the true temperature 
in degrees K) measured by the sonic thermometer. Equation 
(11) is rewritten as: 



T 



SV 



, 2 & 1/2 

''20.067t'' 



( 12 ) 



T is the parameter obtained by measuring t from the son c 

^ V 

thermometer. Being a virtual temperature, this measure- 
ment is subject to humidity effects. 

3 . Discussion of Possible Errors in Wind Velocity 
and Temperature Measurements 

Errors are introduced by humidity, temperature and 
atmospheric pressure which require correcting the readings 
from the sonic anemometer/ thermometer . Detailed discussion 
of errors resulting from misalignment of the array, line- 
averaging, f ini te spatial resolution of the array and noise 
in the detector has been described by Kairaal e_t al^ (1971). 
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As seen in equation (5) and equation (10), the wind 
velocity measured will be a function of temperature, humid- 
ity and atmospheric pressure. From equation (!?.)> however, 
the measured temperature will depend on humidity and atmos- 
pheric pressure only. For hot humid days the effect on 
velocity cannot be ignored. Humidity significantly effects 
the absolute temperature but is negligible as far as fluctua- 
tions in temperature are concerned provided the fluctuations 
in humidity are small. 

The ratio K, between the wind velocity ^^20 

4 

temperature T^q + 293K or 20C and e/p = 0) and the wind 
velocity the temperature is 



K = 



TSV 



V 



20 



SV 

^20 



293 + 6 
293 



SV 



(13) 



where 0^^ = sound virtual temperature measured by the sonic 
thermometer in degrees C. Any wind velocity or temperature 
measured by the sonic anemome t er / t he rmome t e r in an atmos- 
phere o t he r than 20C and zero percent relative humidity 
must be corrected by a factor K (Figure llA, IIB) . 

2 2 

The error introduced by the assumption of C >> V 
is 

Zw = ^ (14) 

C 

and is such that for wind speeds of less than 10 m/sec it 
is zero (Figure 12) . 



Sonic anemome t er / t hermome t er circuitry is set for a base 
temperature of 20C and zero percent relative humidity. 
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V 2 Q = Wind Unit Indica- 
ted Reading 

Vn^ = Wind velocity 

at the temp. 9SV 



dX^Q = Thermometer 
sensitivity 
at 20 degrees C. 

dXg^ = Temperature 
sensitivity 
at temp. T^^ 



Figure llA. Temperature Compensating Constant. . 



te: 



The 0 dial and the meter scale of the TEMP UNIT calibrated in 




Figure IIB. Humidity Compensation Graph. 
Kaijo Denki (1969) 



0 = Terapei ~ture 

= Sound 'irtual 
tempe: iture 
(The temper^iture 
measured by ultrasonic 
thermometer) 

RH = Relative 
humidity 
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4 . Timing Sequence 



A more detailed 
analysis by Souraie and 
Businger (1958) showed 
that velocity and tempera- 
ture fluctuations, within 
certain limits, were pro- 
portional to the fluctua- 
tions in the difference 
and the sum of the time in- 
tervals respectively, A discussion of the method of obtain- 
ing these time intervals is of use in order to understand 
the setting up and calibration of the sonic anemome ter / ther- 
mometer, Only the vertical transmitter-receiver array (that 
used to measure the vertical velocity component and tempera- 
ture) needs to be considered in order to understand the 
triggering and pulse measuring system used in the instrument. 

Pulses are emitted simultaneously from TR__ and TR__ 

w 1 W Z 

at time T, and reach their respective receivers R^^ and R^^ 

after time intervals of t - and t ^ seconds (Figure 13). 

wl w2 VO/ 

The time difference At = t ^ ~ t - is converted 

w w2 wl 

by a wind analog card (within the w unit) and changed into 
a voltage representing the vertical component of the wind 
velocity. The output of R^^ drives the delay card which 
generates the delay time This output synchronizes the 

generation of the next pulse or oscillation at time T’ which 




0 10 20 30 40 SO 



Y C7sec) 



Figure 12. Error by VJind 
Velocity. Kaijo Denki (1969) 



30 



reach their receivers after time intervals t’ ^ and t* ^ 

wl w2 

seconds. Therefore, the pulse spacing is t^^ and t^^ and 
the total time from the emitting of the first pulse, T, to 
the receiving of the second pulse, T’, by is t^^ + t^^^^ 

In terms of the speed of sound in still air, C, and 
the separation, £ , of the transmitter and receiver. 



2 £ 

wl D1 w2 C D1 



(15) 



As previously discussed, the temperature measurement is 
sampled at half the rate of the wind measurement (220 times/ 
sec as compared with 440 times/sec). Therefore, the time, 
t = t^j^ + is related to temperature. Since Atg changes 

very slowly with temperature, a time t ^^2 generated by 



another time delay circuit is subtracted from t + 



The 



time difference, At^ = t + t__ - t_ ^ , is then supplied to a 

U bl DZ 

temperature analog card which changes it into a Atg (the 
mean temperature found on the 0 dial of the temperature 
unit. The dial is operated manually). The time difference 
is also changed into At’g, the fluctuating component, which 
is changed into a voltage that appears on the 0^^ meter 
(Figure 14) . 



The time delay t ^^2 Is adjustable (longer at lower 
temperatures) since it is preferable to keep Atg small. If 
Atg is made too small, however, a sluggish limiting of the 
temperature signal will occur — the effect showing as a 
rounding to flattening of the high temperature peaks. 
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Timing Sequence Schematic. Kaijo Denki (1969) 
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Figure 14. Temperature and Time Relationship. 
Kaijo Denki (1969) 
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Ill, DATA ANALYSIS 



A. GENERAL 

Turbulent processes appear as erratic changes in tem- 
perature and in direction and speed of the wind. A know- 
ledge of the nature of turbulent processes in the atmosphere 
(fluctuations of such parameters as wind velocity, tempera- 
ture and moisture for example) can lead to a better under- 
standing of heat, energy and momentum transfers which occur 
over the ocean. 

The mechanics of turbulence are so complex that it is 
almost traditional now to use statistical methods through 
which an understanding of the mean characteristics can be 
obtained if not in detail. One statistical method which 
has been used with some considerable success is spectral 
analysis. From spectral analysis the distribution of energy 
of a fluctuating, turbulent process (e.g. wind velocity, 
temperature) can be plotted as a function of frequency. 

Analog measurements were recorded of wind velocity and 
temperature fluctuations from the sonic anemometer/ thermo- 
meter, hot wire anemometer and a one mil thermocouple then 
digitized and then spectral analyzed by a Fast Fourier 
Transform method. A general discussion of procedures used 
in the analog-to-digital conversion and spectral analysis 
will be given in the next sections. More detailed discus- 
sion of procedures and theory appear in McKendrick (1971) 
Jones (1971) and Wilson et al (1969). 
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B, KOLMOGOROV’S HYPOTHESIS 



1 . General 

An important consideration in turbulence theory is 
whether the turbulence is anisotropic or isotropic. Aniso- 
tropy exists for large scale eddies and isotropy is thought 
to exist for small scale eddies. Isotropy was described by 
Lumley and Panofsky (1964) as follows: 

”An isotropic field is one in which there is no 
preferred direction. Thus all average functions de- 
scribing the statistics of such a field must remain 
unchanged regardless of a rotation (or reflection) of 
coordinate system. For example, the spectrum of a 
function, (J) (k) , homogeneous in three spatial directions 
which is also isotropic, cannot be a function of the 
separate components of the wave number vector, for a 
rotation to a new coordinate system would change those. 
It can be a function only of the length of the vector, 
for that is the only quantity characterizing ’k’ which 
does remain unchanged under a rotation.” 



turbulence should be isotropic even though it is generated 
from and embedded in anisotropic eddies of low wave number 
If the Reynolds number (ratio of inertial forces to viscou 
forces in the flow) is high enough, turbulence will adjust 



conditions, is attained. Turbulence is adjusted through 
inertial transfer (forces due to fluctuating velocity gradi 
ents) and viscous dissipation. This steady state consists 
of energy transfer by inertial forces into the high wave 
number range that is equal to energy lost through viscous 
dissipation. This predicted behavior of velocity spectra 



Kolmogorov postulated that at high wave numbers 



until a 




steady state, independent of initial 
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was utilized in interpretation and comparison of the sonic 
anemomeme t er , hot wire anemometer and one mil thermocouple 
data. Prediction by Kolmogorov’s hypotheses are discussed 
in more detail in the following paragraphs, 

2 . Kolmogorov’s First Hypothesis: Equilibrium Range 

At high Reynolds numbers ’’small scale” eddies are 
isotropic and in universal statistical equilibrium. Their 
properties depend on dissipation (e) and on kinematic 
viscosity (v) and not directly on time. Motion associated 
with the equilibrium range is uniquely determined by e 
and V , 



Kolmogorov concluded that the energy spectrum 
[E(k)], in units of energy per unit wave number, 



1 s 



E(k) 



= f (^) 



where 



f ('^) ~ undesignated function 







characteristic wave number 



3 , Kolomogor ov ’ s Second Hypothesis: Inertial Subrange 

If the Reynolds number is large enough, the energy 
containing eddy wave number region and the ’’dissipation” 
wave number region should be well separated in wave number 
space. There should exist a range of wave numbers (eddy 
sizes) on the low ”k” side of the dissipation range (but at 
higher wave number than energy containing wave numbers) in 
which the universal statistical steady state holds 
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[r— E(k) = 0] but where there is negligible dissipation of 
o t 

energy, [Figure 15 from Lumley and Panofsky (1964)], This 
range of wave numbers v/here the dominant feature is energy 
flux by inertial forces is called the inertial subrange. 

Since viscosity is not important in this wave number 
range, statistical properties are determined only by dissi- 
pation rate (e) and wave number (k) . Kolmogorov stated that 
energy spectrum should have the following functional rela- 
tion 



E(k) - c 

e = dissipation of turbulent 
kinetic energy 

k = wave number 

C = constant of proportionality 

Since frequency is directly proportional to wave 

number, a log-log plot (base 10) of energy spectra versus 

frequency would yield a curve with a -5/3 slope if the spec- 
tral estimates are within the frequency band of the inertial 
sub range . 



C. ANALOG TO DIGITAL PROCEDURES 
1 . General 

Analog data recorded in the field were digitized 
through use of a Hybrid computer system.^ The Hybrid system 
consists of an analog computer, COMCOR Ci-5000, electron- 
ically connected to a digital computer, XDS-9300. A program 

This system is located in Spanagel Hall, Naval Postgraduate 
School, Monterey, California. 
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Figure 15. Schematic Drawing of Spectral Energy Transfer 
in Free Convection. Lumley and Panofsky (1964) 
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exists for digitizing, using the XDS-9300 and the Ci-5000 
analog computer. An analog patchboard and a logic board, 
which are plugged into the Ci-5000 were also available. 
Digitization was performed in the octal number base onto 
seven track magnetic tapes. 

2 . Considerations Prior to Digitization 

Considerations taken into account prior to digiti- 
zation were noise, D.C. level, analog signal strength and 
total analog data record length, T (corresponding to time). 
All of these were functions of collection technique. 

Noise that contaminates the desired signal is 
usually of high frequency and normally can be low pass fil- 
tered during the data recording procedures. This is the 
preferred method. However, filtering can be done in the 
Hybrid computer facility just prior to digitizing. In this 
study, the analog data were filtered at the Hybrid computer. 
Data from the four sonic anemomet er/ thermometer channels 
were filtered at 31.25 Hz and the hot wire anemometer data 
and one mil thermocouple data were filtered at 100.0 Hz. 
(Filter frequencies should equal one-half of the digital 
sampling rate.) 

The D.C. level was removed from all data during 
digitization. If it were necessary that the D.C. level be 
removed from the data during data recording, a band pass 
filter with a low frequency cut off of less than one hertz 
could have been used. This would preserve some of the low 
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frequencies while eliminating the high frequency noise and 
the D • C . level . 

All analog signals were amplified in order to take 
full advantage of the range of the Hybrid computer system. 
The Hybrid system is able to digitize signals of + 100 volts 
Channels A and B of the sonic anemometer/ thermometer and 
the hot wire anemometer were amplified by 50; Channels W 
and T of the sonic anemometer/ thermometer and the thermo- 
couple were amplified by 100. Amplification factors were 
taken into account, in this study, in the Fast Fourier Trans 
form program. 

Consideration had to be given to obtaining an opti- 
mum analog data record length, T, during the field experi- 
ments. In optimizing this total record length for the 
frequencies of interest, the relationship between T and the 
shape of the convolved wave forms (Fourier Transform) must 
be considered. Jones (1971) described fully the effects of 
record lengths on the shape of the convolved wave forms. 

In general there is some optimum record length where the 
difference between the convolved spectrum and true spectrum 
is very small. Although Jones (1971) did not find an opti- 
mum record length, he found that increasing a record from 
140 seconds to 360 seconds had no significant effect on 
spectral shape between 10 Hz and 1000 Hz. Data collected in 
this study had a record length of 30 minutes. 

Low frequency resolution is also tied directly to 
the length of the analog record. According to Rayleigh^s 
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Criterion, it is hard to distinguish low frequencies (f) in 
the interval 0 f ^ 1/2T (where T = analog record length 
in seconds). When dealing with a finite record the original 
signal is truncated, A long record narrows the frequency 
resolution. Record lengths of 30 minutes gave a low fre- 
quency resolution of approximately 0,0003 Hz although sta- 
tistical reliability would not be expected until 0,003 Hz, 
Another factor which had to be considered prior to 
digitizing was that of aliasing. Aliasing is dependent on 
the folding (Nyquist) frequency which is 1/2 the sampling 
frequency. In order to prevent the”folding back” of energy 
from higher frequencies into the lower frequencies of 
interest, a sampling rate greater than twice the highest 
frequency desired must be chosen. The folded back energy 
will be centered around the Nyquist frequency. Sonic anemo- 
meter/thermometer data from Experiment I was run at two 
different sampling rates. Sampling rates of 390.5 Hz and 
62,5 Hz gave aliasing frequencies of 195.25 Hz and 31.25 Hr 
respectively. The hot wire anemometer data and the thermo- 
couple data from Experiment II were sampled at 200.0 Hz which 
gave an aliasing frequency of 100.0 Hz, Experiment II sonic 
anemometer/ thermometer data were sampled at 62,5 Hz only. 

D, CONVERSION OF SEVEN TRACK TO NINE TRACK 

Digitization of analog data on the above system produced 
a seven track tape in the octal number case. To utilize 
this data in the IBM-360/67 computer system, digitized data 
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had to be converted to a nine track tape in the hexadecimal 
number base. This conversion process was done using an 
available program, CONVERT. Complete listings and Job Con- 
trol Language (JCL) card setups for the CONVER.T program can 
be found in theses by McKendrick (1971) and Jones (1971), 
and in a report by Raney (1971). 

E. FOURIER COEFFICIENT COMPUTATIONS 

Fourier coefficients were obtained by using the FTOR 
program, one of two UBC time series programs currently on 
file in the IBM-360/67 subroutine library at the Naval Post- 
graduate School. 

Transformation of digitized signal data from the time 
domain into the frequency domain was done using a finite 
form of the Fourier Transform. Data were read off the con- 
verted nine track tape, a record at a time. The Fourier 
coefficients were computed and were then stored on a sepa- 
rate tape. 

The longer the. record, the closer together in frequency 
will be the spectral estimates because the value of the 
frequency is given by (n/T) . That is, the discrete frequency 
bandwidth is inversely proportional to the record length. 

Due to the prohibitive amount of computer time used to 
compute coefficients, the direct finite Fourier Transform 
of a series of data points was not used. The Fast Fourier 
Transform of Cooley and Tukey (1965) was used and it is a 
variation of the direct method. The Fast Fourier Transform 
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is incorporated into the FTOR program resulting in consider- 
ably less computer time for Fourier coefficient computation. 
This method yields an improvement in excess of five orders 
of magnitude in computational time (McKendrick 1971) over 
the direct method, 

13 

The FTOR program transforms a maximum of 2 (8192) data 

points. This maximum is set by dimension limitations in the 

program. Turbulence data normally includes many more data 

points than this so the data points are read in as records 

consisting of integral powers of two (up to a maximum of 
l3 

2 ) and these records are transformed consecutively, 

F, UBCSCOR PROGRAM 



Once the Fourier coefficients have been computed, the 
second UBC program, SCOR, converts them to power spectral 
densities . 

Spectral values E., (n) are computed from the complex 
Fourier coefficients by forming the product of the complex 
coefficient [Aj (n) ' = a^ (n) - i b^ (n) ] with a conjugate 

(n) = aj^(n) + i bj^(n)] where j and k could be the same 
or different. 

Spectral values were obtained by the following computa- 
tions 



Spectrum density: j = k. 



E . . (n) 
22 



A. (n) ‘A. (n) 
= 2 2 
2 An 



where An = ^ 




+ 



n 



T 

2 
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Cross spectrum ; j ^ k, 

A (n) * A, (n) 

2An" ' 

■ f 

The real and imaginary parts of the cross spectrum are 
the cospectrum and quadrature spectrum respectively and 
were denoted as; Cospectrum : 

E (n) = (a c + b d ) ^ 

3k n n n n 2 

Quadrature spectrum : 



E (n) =(bc -ad)— 
jk nn nn2 



The spectra and cospectra are computed for each fre- 
quency in successive records and then these records are 
averaged consecutively. The program also divides the spec- 
tral densities into 32 bandwidths and averages them across 
the bandwidths as well as along the records to obtain a 
mean and a variance of the quantities. 

In the SCOR program one can specify the number of 
channels to be analyzed, the number of records in a file, 
the start point for analysis, selection of spectra and/or 
cospectra, bandwidth desired (logarithmic or linear), the 
type of graphic plot desired and for this study, machine-set 
coordinates for the graphic output. Programmer-set-co- 
ordinates also can be used. The program computes, the 
spectra and/or cospectra for each record in sequence until 
all records requested have been analyzed or the ”end-of- 
file” mark is reached on the tape. 
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IV. FIELD EXPERIMENTS 



A. PRELIMINARIES 

The sonic anemometer/ thermometer used in the field ex- 
periments had been in storage for a period of two years. 

The first consideration was to thoroughly check the sonic 
anemome t er / t he rmome t er electronically. It was checked out 
in accordance with manufacturer’s instructions [Kaijo Denki 
(1969) , Appendix D and E which included voltage and wave form 
checks at certain critical points]. A large number of tran- 
sistors in both the power supply and junction box had to 
be replaced . 

Once the sonic anemometer / thermometer was checked elec- 
tronically, wind tunnel tests were performed. Simple, 
operational tests in a controlled environment were made to 
insure that the sonic anemometer /thermometer was measuring 
wind speed and temperature correctly. Once satisfactory wind 
tunnel measurements were made, field experiments were 
conducted . 

B. OBJECTIVE OF EXPERIMENTS 

The objective of the field experiments was to gather com- 
parative turbulence data from the sonic anemometer /thermo- 
meter, the hot wire anemometer and a one mil thermocouple. 

All three sets of data were then similarly spectrally 
analyzed. 
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EXPERIMENTS 



C . 

1 . Installation of Instruments 

The procedures described in this section applies to 
both experiments. All equipment was installed at the Naval 
Postgraduate School beach lab , in di ca t ed in Figure 16 and 
Figure 17. The sonic probe, hot wire probe and one mil ther- 
mocouple were mounted on a pole at a height of three meters 
(Figure 18) and positioned so that all sensors were facing 
directly into the wind flow which was blowing off the water. 
During the experiments, all probes were kept into the wind 
such that the velocity measured on channels A and B were 
equal. A hand held wind vane was used to check on the wind 
direction in order to maintain proper probe orientation. 

All recording devices , power supplies and amplifiers 
were installed inside the beach lab (Figure 19). Analog 
data were recorded on a SANGAMO FM tape recorder model 3562 
at a speed of 1-7/8 inches per second. This recorder uses 
one inch magnetic tape and has 14 tracks available for data 
and two edge tracks available for voice recording. 

The hot wire anemometer used for velocity compari- 
sons was a TSI model 1050. A copper constantan, unsheathed, 
fine gauge Omego Engineering Company ( S ST S SP - 0 2 OE- 6 ) thermo- 
couple was used for temperature measurements. 

2 . Experiment 1-14 December 1972 

Wind speed and temperature measurements were made 
during the late afternoon on 14 December 1972. Visibility 
was unlimited. Winds were from the west to northwest, off 
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Figure 16 . 



Arrangement of Probes for Experiments 



I and II . 
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F igur e 17 . 



Closeup 



of Probe Arrangement. 
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Figure 18 . Probe Stand Orientation. 
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Figure 19. Inside Installation of Equipment. 
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the ocean, at an average speed of 2-4 ra/sec. Mean air tem- 
perature for the period (1430-1700 PST) was 52°F. 

Three separate 30 minute analog data records were 
recorded (1500 to 1530, 1546 to 1616 and 1620 to 1650 PST). 

Data on all channels were recorded without amplification. 
During the measurements, spikes were observed on all channels. 
The spikes were found, after this experiment, to have been 
caused by faulty transistors in the junction box. Examples 
of these spikes appear in the strip chart recording in 
Figure 20. 

Hot wire anemometer and thermocouple data were not 
recorded because the record amplifier on the hot wire channel 
of the tape recorder failed and the thermocouple malfunc- 
tioned. 

The three separate analog data records from the 
sonic anemometer/ thermometer data were digitized and spectral 
analyzed . 

3 . Experiment II - 27 February 1973 

This experiment was similar to experiment I except 
that hot wire anemometer data and thermocouple data were 
successfully recorded. Channel B of the sonic anemometer/ 
thermometer was not operating because of a short in the con- 
nector cable on that channel. This malfunction did not 
affect the other channels. 

Again, wind speed and temperature data were obtained 
during three separate 30 minute periods (1030 to 1100, 1345 
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Figure 20. 



Strip Chart Recording of the Sonic 
Channels Indicating Spikes. 



and W 
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to 1415 and 1530 to 1600 PST), It was an overcast sky. 
Visibility was approximately six miles. Winds were from the 
west, off the ocean, with an average speed of 4-8 m/sec. 

Mean air temperature ranged from 58°F in the late morning 
period to 60®F in the early afternoon period, 

During the recording periods, all but one of the 
data channels were amplified using Preston 8300-XWB ampli- 
fiers. The amplification factors were two for the sonic 
W and T channels, one for the sonic A channel, twenty for 
the hot wire channel and one thousand for the thermocouple 
channel. The hot wire data were filtered before recording 
using a Krohn-Hite model 3750 bandpass filter low set at 
0.2 Hz and high set at 200,0 Hz. The three separate 30 
minute analog data records were digitized and spectrally 
analyzed. 
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V, 



RESULTS AND CONCLUSIONS 



A. GENERAL PERFORMANCE OF SONIC ANEMOMETER/ THERMOMETER 
In general the sonic anemometer/ thermometer satisfac- 
torily measured the wind and temperature fluctuations during 
the field experiments, but it did not operate at its peak 
potential because of various electronic problems. For 
example, spiking occurred on all channels during both ex- 
periments I and II and a short occurred in channel B during 
experiment II. These could have been a result of the instru- 
ment being idle for two years. 

Spiking proved to be a continuous problem during the 
field experiments even though numerous transistors in the 
system had been replaced. 

B. SPECTRA PLOTS 
1 . General 

Sonic anemometer/ thermometer data recorded on 14 
December 1972 (experiment I) were digitized at a sampling 
rate of 390.5 Hz (aliasing frequency 192.25 Hz) and also 
at 62.5 Hz (aliasing frequency 31.25 Hz. 

No editing or filtering was performed on any of the 
channels during digitization at the higher sampling rate. 

The higher sampling rate (390.5 Hz) enabled the examination 
of a frequency band from 0 Hz to approximately 100 Hz. 

The analog data tape was edited prior to digitizing 
at the lower sampling rate. This led to the elimination of 
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the first time period which had the heaviest spiking. The 
remaining two time periods were digitized at 62,5 Hz and 
filtered at 31.25 Hz. The lower sampling rate shifted the 
frequency band examined (approximately 0.5 Hz to 40 Hz) 
towards the lower end of the frequency spectrum. The spec- 
tral results indicate both the energy containing frequency 
band as well as the lower frequency end of the inertial sub- 
range. 

Hot wire data were not taken during experiment I 
because of a recording amplifier failure. Two samples of 
hot wire data taken on a day that had similar synoptic 
characteristics (air temperature and wind) as the day of 
experiment I are included for comparison purposes. They 
are considered representative samples of data which would 
have been taken by a hot wire anemometer system on that day. 
These results represent a frequency band of from 0.01 Hz 
to approximately 100 Hz. The data for both of these results 
were digitized at a rate of 200.0 Hz (aliasing frequency 

100.0 Hz) and filtered at 100.0 Hz. 

On 28 February 1973 (experiment II), sonic anemo- 
meter/thermometer, hot wire anemometer and one mil thermo- 
couple data were collected during the three 30 minute time 
periods. The sonic data were digitized at a sampling rate 
of 62.5 Hz and filtered at 32.25 Hz. The hot wire and 
thermocouple data were digitized at a sampling rate of 

200.0 Hz and filtered at 100.0 Hz. Sonic channel B was 
inoperat ive . 
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2 • Sonic Anemome t e r /The rmome t er Data Sampled at 

390.5 Hz - 14 December 1972 (Experiment I) 

Figure 21 to Figure 32 are for the three 30 minute 
periods from 1500 to 1700. A -5/3 slope exists on the log- 
log plot in the frequency band from 1.0 Hz to 10.0 Hz for 
all channels. Spectra from channels A, B, W have -5/3 
slopes which continue beyond 10.0 Hz but the T spectra 
slopes drop off more rapidly than -5/3 before 10.0 Hz. 

3 . Sonic Anemometer /Thermometer Data Sampled at 

62.5 Hz - 14 December 1972 (Experiment I) 

Figure 33 to Figure 40 are from the two edited 
periods of data which were filtered and then sampled at the 
lower rate. Spectra from channels A and B in both periods 
have definite -5/3 slopes within the 1.0 Hz to 10.0 Hz band 
with some departure around 8.0 Hz. For both periods, W 
channel spectra have -5/3 slopes up to 6.0 Hz with rapid 
deterioration beyond. For the second period, the T spectrum 
has a drop off from a -5/3 slope at approximately 4.0 Hz but 
for the third period, the T spectra has a definite -5/3 
slope throughout the 1.0 Hz to 10.0 Hz frequency band. 

4 . Hot Wire Anemometer Data Sampled at 200.0 Hz 

- 10 January 1973 

Figure 41 and Figure 42 are from hot wire anemometer 
data taken on 10 January 1973. Two periods are discussed 
and each period was approximately 50 minutes in length. A 
distinct -5/3 slope is found in the frequency band from 
1.0 Hz to approximately 40.0 Hz for both periods of data. 
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Figure 21, Sonic Channel Data Observed 14 December 

1972 1st Period (Experiment I) High Sample Rate, 
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Figure 22. Sonic "B" Channel Data Observed 14 Decenber 
1972 1st Period (Experiment I) High Sample Rate. 
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Figure 23. Sonic Channel Data Observed 14 December 

1972 1st Period (Experiment I) High Sample Rate. 
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Figure 24. Sonic Channel Data Observed 14 December 

1972 1st Period (Experiment I) High Sample Rate. 
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Figure 25 . Sonic "A" Channel Data Observed 14 December 
1972 2nd Period (Experiment I) High Sample Rate. 
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Figure 26 , Sonic "B" Channel Data Observed 14 December 
1972 2nd Period (Experiment I) High Sample Rate. 
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Figure 27. Sonic Channel Data Observed 14 December 

1972 2nd Period (Experiment I) High Sample Rate. 
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Figure 28. Sonic "T" Channel Data Observed lA December 
1972 2nd Period (Experiment I) High Sample Rate. 
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Figure 29. Sonic Channel Data Observed 14 December 

1972 3rd Period (Experiment I) High Sample Rate. 
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Figure 30. Sonic "B" Channel Data Observed 14 December 
1972 3rd Period (Experiment I) High Sample Rate. 
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Figure 31. Sonic "W" Channel Data Observed 14 December 
1972 3rd Period (Expeiiment I) High Sample Rate. 
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Figure 32, Sonic "x" Channel Data Observed 14 Deceaiber 
1972 3rd Period (Experiment I) High Sample Rate. 
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Figure 33. Sonic Channel Data 

1972 2nd Period (Experiment I) Low 
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Figure 34. Sonic Channel Data Observed 14 December 

1972 2nd Period (Experiment I) Low Sample Rate. 
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Figure 35. Sonic W" Channel Data Observed 14 December 
1972 2nd Period (Experiment I) Low Sample Rate. 
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Figure 36. Sonic **T” Channel Data Observed 14 December 
1972 2nd Period (Experiment I) Low Sample Rate. 
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Figure 37 . Sonic "A" Channel Data Observed 14 December 
1972 3rd Period (Experiment I) Low Sample Rate. 
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Figure 38. Sonic "B" Channel Data Observed 14 Decenber 
1972 3rd Period (Experiment I) Low Sample Rate. 
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Figure 39. Sonic **W” Channel Data Observed lA December 
1972 3rd Period (Experiment I) Low Sample Rate. 
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Figure 40. Sonic "x" Channel Data Observed 14 December 
1972 3rd Period (Experiment I) Low Sample Rate. 
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Figure 41. Hot Wire Data Observed 10 January 1973 
1st Period . 
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Figure 42, Hot Wire Data Observed 10 January 1973 
2nd Period . 
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5 . Sonic Anemometer/Thermometer Data Sampled at 
62.5 Hz - 28 February 1973 (Experiment II) 

Figure 43 to Figure 51 are spectra from the sonic 

anemometer / thermometer data recorded on 28 February 1973. 

Channel B was inoperative during the period. Spectra from 

channels A and W show a definite -5/3 slope in the 1.0 Hz 

to 10.0 Hz frequency band. A -5/3 slope in a very narrow 

band of frequencies of approximately 1.0 Hz to 4.0 Hz 

appears in the T spectra. 

6 . Hot Wire Anemometer Data Sampled at 200.0 Hz 
- 28 February 1973 (Experiment II) 

Figure 52 and Figure 53 are hot wire spectra from 
the second and third periods recorded on 28 February 1973. 
The first period of data was lost due to a faulty magnetic 
tape. Spectra from both periods have definite -5/3 slopes 
in the region 1.0 Hz to approximately 91.0 Hz. 

7 . Thermocouple Data Sampled at 200.0 Hz - 
28 February 1973 (Experiment II) 

Figure 54 and Figure 55 are thermocouple spectra 
from the second and third periods recorded on 28 February 
1973. The first period of data in this case was also lost 
due to the faulty magnetic tape. From the second period 
plot the -5/3 slope starts at approximately 4.4 Hz and goes 
just beyond 10.0 Hz. From the third period plot, a 
definite -5/3 slope exists in the range 1.0 Hz to 10.0 Hz. 

C . CONCLUSIONS 

The sonic anemometer/ thermometer system proved to be a 
reliable, rugged and easily installed instrument for the 
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Figure 43. Sonic "A" Channel Data Observed 28 February 
1973 1st Period (Experiment II). 
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Figure 44. Sonic "W" Channel Data Observed 28 February 
1973 1st Period (Experiment II). 



81 



lOGLO Sr’tGTRUM ( /HGRFZj 




Figure 45. Sonic "x" Channel Data Observed 28 February 
1973 1st Period (Experiment II). 
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Figure 46, Sonic Channel Data Observed 28 February 

1973 2nd Period (Experiment II). 
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Figure 47. Sonic "W" Channel Data Observed 28 February 
1973 2nd Period (Experiment II). 
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Figure 48. Sonic "T" Channel Data Observed 28 February 
1973 2nd Period (Experiment II). 
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Figure 49. Sonic "A** Channel Data Observed 28 February 
1973 3rd Period (Experiment II). 
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Figure 50, Sonic ”W" Channel Data Observed 28 February 
1973 3rd Period (Experiment II). 
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Figure 51. Sonic "T" Channel Data Observed 28 February 
1973 3rd Period (Experiment II). 



88 



) Untii'J3dS OTGO 



u 




Figure 52. Hot Wire Data Observed 28 February 1973 
2nd Period (Experiment II). 
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Figure 53. Hot Wire Data Observed 28 February 1973 
3rd Period (Experiment II). 
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Figure 54. Thermocouple Data Observed 28 February 1973 
2nd Period (Experiment II). 
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Figure 55. Thermocouple Data Observed 28 February 1973 
3rd Period (Experiment II). 



92 



measurement of the three components of the wind and also air 
temperature. Its frequency response was found to be from 
0.0 Hz to approximately 12.0 Hz. 

The sonic anemometer/ thermometer measurements did 
not appear to be affected by weather elements and because it 
uses the speed of sound in air it has inherent calibration. 
The response of the sonic anemometer/ thermometer is linear 
and so is its sensitivity. It is thus able to be used to 
describe vertical and horizontal transfers of energy, heat 
and momentum. 

The ~5/3 slope of Kolmogorov’s inertial subrange 
was easily recognizable on spectral plots of turbulence data 
measured by the sonic anemometer/ thermometer . Thus the com- 
putation of the rate of dissipation of kinetic energy con- 
stant can be made which in turn means that heat and momentum 
transfer computations can be made. 
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APPENDIX A 



Definition of Terms Used in Analog-Digital Data Processing 

The terms used in the analog-to-digital process are de- 
fined here based on current usage in the Hybrid computer 
center and based on that utilized by the three time series 
programs for the spectral analysis of digitized data. There 
are some contradictions of terms between McKendrick (1971) 
and Jones (1971) and between what is currently used. 

The word "record” is often used in both analog and digi- 
tal discussions as meaning a finite length of data in 
seconds. In the analog sense, "record" means the entire 
length of data which is to be digitized. In the digital 
sense, "record" means a short length of data based on the 
sampling rate and the number of samples per record selected 
by the programmer. In McKendrick (1971) and Jones (1971), 
digital "record" is referred to as "block". The word 
"block" is not used in this thesis as "record" is what is 
utilized in all three programs (although the print outs use 
blocks). A digital "record" is limited in the number of 
samples it can contain because of Hybrid computer memory 
limitations. This limit is that the number of samples t imes 
the number of channels simultaneously digitized must be 

equal to or less than 6144 channel -samples . As an example, 

record 

consider the following hypothetical problem: 
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Analog record length « 3,000 secs 

Number of analog channels containing data = 5 channels 

Sampling rate = 100 s amp le s 

sec 

Number of channels digitized 

simultaneously = 5 channels 

Samples per record selected = 1024 samp les 

record 

1024 — ~ X 5 channels = 5120 channel-saraole ^ 

record ^ ‘ < 6144 

record — 

samples 1 sec seconds 

record 100 samples record 

Digital record length = 10.24 seconds 

record 

Note that the value for samples for record was chosen as 

an integral power of two , i.e. 2^^, The reason is that the 

max imum number of data points (samples) that can be stored 

in the IBM 360/67 computer memory under the FTOR program are 
13 

8192 or 2 samples. Therefore, for efficiency purposes t ae 
FTOR program was designed to read in samples per record wl .ch 
contain integral powers of two. 

One other commonly used word is "file.” A ^'file" is a 
finite length of digital data containing any number of 
records from one on up. The length of a "file" is predeter- 
mined during the digitization phase and is based on the 
number of records desired per "file." An "end of file" 
mark is placed on the seven track tape after the desired 
number of records have been computed. 
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